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Abstract 

We calculate the dynamics of a nanomechanical oscillator (NMO) coupled capacitively to a 
Cooper pair box (CPB), by solving a stochastic Schrodinger equation with two Lindblad operators. 
Both the NMO and the CPB are assumed dissipative, and the coupling is treated within the rotating 
wave approximation. We show numerically that, if the CPB decay time is smaller than the NMO 
decay time, the coupled NMO will lose energy faster, and the coupled CPB more slowly, than do 
the uncoupled NMO and CPB. The results show that the efficiency of energy loss by an NMO can 
be substantially increased if the NMO is coupled to a CPB. 



PACS numbers: 



I. INTRODUCTION 



Observing quantum effects in a mechanical object has been a long-sought goal, because 
such effects would represent a macroscopic manifestation of quantum mechanics[l|. Nanome- 
chanical oscillators (NMO's) with frequencies uj/{2ti) as high as 1 GHz have recently been 
realized in the laboratory!^. It is possible to observe quantum effects in such objects, pro- 
vided that they are cooled to temperatures T such that ksT is near or slightly above ^jsl- 
For a realistic NMO frequency of 500 MHz, this corresponds to a T of a few mK. 

In this paper, we numerically analyze energy loss by a NMO capacitively coupled to a 
dissipative Cooper pair box (CPB). A CPB consists of a small Josephson junction whose 
Josephson energy is much smaller than its charging energy, and which is voltage-biased 
so that only the two lowest energy levels are experimentally accessible. Naik et a/. [4] have 
recently suggested that an NMO could be cooled by capacitively coupling it to a CPB. In the 
present work, we numerically demonstrate an increased rate of energy loss by this mechanism, 
by solving a time- dependent Schrodinger equation with a stochastic term. We also relate 
the rate of energy loss to the coupling constant and the NMO and CPB dampiiig rates. 
Several studies of the coupled NMO/CPB system have already been carried out js], 0,0], as 
have studies of the formally similar system consisting of a CPB couplied to a single-mode 
electromagnetic cavity j^. If both the CPB and the cavity are dissipative, this system is 
closely analogous to the one studied here ^ . 

II. FORMALISM 

Dissipation in quantum systems was treated by Caldeira and Leggett, who studied the 
suppression of quantum tunneling out of a metastable state by dissipation [Oj. Later, Leggett 
et al. studied the dynamics of a dissipative TLS [l^ . In both of these approaches, the dissi- 
pation is treated as arising from coupling to a bath of harmonic oscillators. An alternative 
method starts from Lindblad's quantum-mechanical master equation for the density opera- 



tor of an open system 



111] . An equivalent approach is to start from a pure-state Schrodinger 



equation, which contains a stochastic term to represent the effect of dissipation Here, 
we use this approach to treat energy loss by an NMO coupled to a CPB. 
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We assume that the Hamiltonian for the coupled NMO / CPB system is given by 



H = HcpB + Hnmo + Hj + + Hj, (1) 
HcPB = ^Eciug - ^)a, - ^Eja^, (2) 

Hnmo = hwo{a^a + ^), (3) 
Hi = hg{a) + a)az. (4) 

Here az and are the Pauh spin matrices, which act on the two charge states of the CPB; 
Hg = {CbVb + CgVg)/2e, where Cb and Vb are the bias capacitance and voltage of the CPB; Cg 
and Vg are the capacitance and voltage between the NMO and the CPB; Ec and Ej are the 
Coulomb and Josephson energies for the CPB; ujq is the frequency of the fundamental mode 
of the NMO; a and are the usual annihilation and creation operator for the fundamental 
mode of the NMO, and g = —AEc{CgVg)/S.Xzp/{2ehd) represents the coupling constant for 
capacitive coupling between the NMO and the CPB. The zero-point displacement of the 
NMO is given by Ax^p = where m is the effective mass of the NMO and d is 

the distance between the NMO and the CPb[i3|. The last two terms in H represent the 
effects of dissipation, if^ represents the dissipation in the NMO, which is assumed to be 
characterized by a decay rate k = cuq/Q, where Q is the quality factor of the NMO. H^ 
represents the dissipation in the CPB, which is described by a decay rate 7. We also assume 
that only the fundamental mode of the NMO is coupled to CPB. The precise forms of H,^ 
and H^ are given below. 



The Hamiltonian (1), in the absence of dissipation, has been discussed, e. g., in Refs. 



(which shows a sketch of a typical experimental configuration), [6] and [15|. We assume that 
the CPB is biased to = 1/2, where it would be degenerate in the absence of Josephson 
coupling, and work in the coordinate system obtained from that of eq. (jl]) by a rotation 
of 7r/2 about y-axis. In this coordinate system Hcpb = {Ej/2)az and Hj = hg{a + a'^)ax- 
Moreover, in this coordinate system, the two eigenstates of Hcpb represent symmetric and 
antisymmetric superpositions of charge states. In addition, we will make the rotating wave 
approximation (RWA) to obtain the Hamiltonian used in all our calculations. This Hamil- 
tonian is given by 

H = ~^'^z + ^0 ^ct^a + 2^ ~'~ %('^^^^ + a^a) + H^ + if^, (5) 
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where Q = Ej/h and = {a^ ± iay)/2. Thus, we neglect the counterrotating terms 
h£(a^cr~^ + aa~) If ~ and the couphng is sufficiently weak, the RWA is believed justified 
131]. For this reason, we believe that the RWA is also reasonable when the two damping 
terms are included. 

We include dissipation in both the CPB and the NMO by using a stochastic time- 
dependent Schrodinger equation. This approach originates in the von Neumann equation 
for the total density matrix p of any system: p = —j^[H,p]. The Hamiltonian of eq. ([5]) 
consists of two parts: the CPB plus NMO, described by the first three terms, and the dissi- 
pative parts, described by iJ^ and H^, which represent the damping. If the evolution of p is 
considered to be Markovian, and if the coupling between the open system and the baths is 
weak, the reduced density matrix p.o,,,^ of the CPB plus NMO can be shown to be described 



by the following master equation 



16|: 



Psys — ~t[H — — H^, Psys] — - ^ {Ll^Lmpsys + PsysLl^Lm — 2LmpL\^)- (6) 



m={K,7} 

Here the Lm are so-called Lindblad operators, which are in general non-Hermitian and which 
describe dissipation in CPB plus NMO system. In the present work, we include only the 
two Lindblad operators = ^/na and = ^/^a" , which represent dissipation within the 
NMO and the CPB. The inclusion of only these two Lindblad operators is appropriate if T 
is sufficiently low (typically /c^T <^ huo), and if other dissipative processes can be neglected. 

Eq. ([6]) is a matrix equation for the density operator p^ys- If we truncate the state 
space by including only the lowest vibrational states, then p^ys is a square matrix of size 
2N X 2N. Eq. ([6]) then leads to (2N)^ coupled differential equations for the elements of p^y.^ . 



An alternative computational scheme is to use a stochastic pure state representation 12 



li 



for the master equation. In this representation, eq. (E]) reduces to a stochastic equation of 



motion for the state vector given by 
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m 

- ^ E {^^rnLm + (^L)^ (Lm)^ ' 2 (L^) ^ L^) |^) dt, (7) 
m 

where {Lm)^ is the expectation value of Lm in state and Idip) is the change in in 
time dt. In eq. ([7]), the first sum represents random fluctuations due to the interaction of 
the system with the baths, while the second denotes the (non-random) drift of the state 



vector due to those baths. The d^m are independent complex differential random variables 
representing a complex normalized Wiener process, whose ensemble averages satisfy 



n 



(8) 
(9) 



where the overbar denotes an ensemble average over realizations of the random variables 



Eqs. ([7]) represent only 2N coupled differential equations for the components of rather 
than {2NY as in the density matrix formulation. The price paid is that the differential 
equations are stochastic, and thus must be averaged over many realizations. But such 
stochastic equations can be solved very efficiently, even including this averaging. Thus, we 
use this stochastic approach in the following. 

To see how this procedure can lead to faster energy loss by the NMO, suppose that the 
initial state vector is \ip) = \<y,n), where a{= or 1) indicates an eigenstate of Hcpb and 
n represents the initial number of excitations in the NMO. Thus, the NMO is assumed to 
be initially in a Fock (or number) state. As time progresses the NMO and CPB become 
entangled, and the wave function spreads out over an increasing part of the accessible 
2N-dimensional Hilbert space. After some time t, \'4'{t)) is a superposition of many states 
|a,m), with m < n + 1; the weight of each component is determined by eq. ([7]). Several 



III. NUMERICAL RESULTS 

We have solved eq. ((71) at resonance for different values of k and 7, using the value 
g = O.ltuo, by straightforward Euler integration, using a time step dt = 5 x 10^^ns= 2.5 x 
10^^(27r/co'o). From these solutions, we compute various averages of interest, such as {(Tzit)) 
and {a^a{t)), as functions of time t. The triangular brackets and overbars denote quantum- 
mechanical and noise averages, respectively. For conciseness, we denote these quantities 
simply a^{t) and n{t). The results are averaged over 100 realizations of the noise, which 
appear sufficient to give reasonably smooth results. We carry out calculations assuming 
f)/(27r) = cuo/(27r) = 500 MHz, g/{2n) = -50 MHz, n = 0.0005 (ns)-iand 7 = 0.05 (ns)-^ 
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These values of k and 7 correspond to decay times of the uncoupled NMO and CPB of 
2000 ns (corresponding to a Q factor of 10'^) and 20 ns, respectively, as used in Ref. jl8 |. 
A much larger decay time for the CPB has been achieved in some more recent experiments 
(see, e. g., Wallraff et a/, [l^, who find an energy relaxation rate of around Tyusec). However, 
one does not want such a long decay time for the CPB, since this will reduce the rate at 
which the NMO loses energy when coupled to the CPB. We have, therefore, carried out our 
calculations using the original parameters of Ref. jl8 |. 

We first carried out calculations of the uncoupled CPB and NMO, using the stochastic 
Schrodinger equation, starting from the state |1,20). In this case, the CPB and NMO each 
decay exponentially with decay rates 20 (ns)~^ and 2000 (ns)~^, respectively, as shown in 
Figs. 1 and 2. Next, we turned on the coupling g, so that g/{2'K) = —50 MHz, with results 
shown in Figs. 1 and 2 for o^{t) and n{t). 



IV. DISCUSSION 



We now discuss the results. In both cases, we use 7 ^ k of the NMO, as in experiment 



0, 



18j. Fig. 1 shows that, in the presence of coupling, the decay rate of the CPB is substantially 



reduced, relative to the uncoupled CPB. The reduction in the decay rate is greatest when 
the damping in the NMO is smallest. Besides this reduction in decay rate. Fig. 1 shows the 
expected Rabi oscillations in for all three calculations involving nonzero g. These Rabi 
oscillations, of frequency uj^i oc \/n + lq[2o|. occur between the states and |0,n + 1). 

They are the well-known results of solving the time- dependent Schrodinger equation for the 
present Hamiltonian in the absence of dissipation; our calculations show that they still occur 
with dissipation present. 

By contrast. Fig. 2 shows that, when the coupling is turned on, the decay rate of the 
NMO is mcreasec? relative to that of the uncoupled NMO. Furthermore, the decay rate is the 
largest when the damping rate of the CPB is largest. Once again, we see Rabi oscillations 
in n(t). 

The results of both Figs. 1 and 2 can be understood qualitatively when 7 ^ k. The 
coupling allows energy to be transferred periodically between the CPB and the NMO. Since 
the NMO has a much lower decay rate than that of the CPB, this new channel should 
reduce the effective decay rate of the CPB, as we observed numerically. Similarly, this 
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transfer should increase the effective damping rate of the NMO, again as seen numerically. 
Furthermore as the damping rate 7 of the CPB increases, the effective damping rate of the 
NMO should also increase, again as seen numerically. 

Finally, in Fig. 3, we show how this behavior depends on the assumed initial Fock state 
of the NMO. In the three cases shown, we start from state |l,n) with four different initial 
Fock states, n(0), of the NMO. n{t)/n{0) decays roughly exponentially with time, with 
superimposed Rabi oscillations. The exponential decay rate is consideably smaller for n{0) = 
1 than for the other cases. This difference occurs, we believe, because, in the absence of 



damping, the coupling between the CPB and NMO is proportional hg^/n + 1 (this is the 
splitting between states \l,n) and |0,n+ 1) in the RWA in the absence of damping). 

We now discuss the possible relation between the present results and energy loss by a real 
NMO. We have considered the time-dependent NMO/CPB system, starting from the state 
|l,n). Our approach permits us to calculate the time needed for the NMO, starting from 
this excited state, to approach its equilibrium temperature, given that the intrinsic damping 
within the NMO is very small. If the NMO is coupled to a much more heavily damped 
CPB, we find that the equilibration process is greatly speeded up. This result might be 
useful in designing ways of rapidly equilibrating an NMO, initially in an excited state, to a 
low ambient temperature. 

We also find that the rate of equilibration of the CPB is reduced when it is coupled to 
an NMO. Our results are consistent with those of Trees et a/. {21}], who used a very different 
approach. These workers studied a current-biased Josephson junction capacitively coupled 
to an oscillator, in the presence of dissipation, using the formalism of Caldeira and Leggett ^ , 
and found that the CPB is less damped when it is coupled to an NMO. Our results show 
the same behavior, using a formally quite different stochastic differential equation. We also 



find an additional result, not shown in Ref. 



2l| : for a given coupling constant, the CPB is 



damped most slowly when the damping of the NMO is smallest (see Fig. 1). 

It should also be pointed out that we have included only two Lindblad operators in our 
calculation, namely, those involving the constants k and 7. Other Lindblad operators could 
readily be included, such as ones describing finite temperature and pure dephasingjs^. In the 
present calculation, k and 7 parametrize energy loss from the NMO and the CPB to their 
thermal environments. Another point is that the regime we have considered, in which the 
CPB is more heavily damped than the NMO, is not necessarily the regime in which a CPB is 
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typically studied. When one wishes to use the CPB as an element in a quantum computing 
geometry, it is usually desirable to have a CPB with as httle dissipation, and as long a 
decoherence time, as possible. In the present work, by contrast, we are envisioning a setup 
where the CPB is being used as a means of rapidly extracting energy from a high-Q NMO, 
so as to observe the expected quantum effects in the NMO. As shown by our calculations, 
this requires a CPB which is more heavily damped than the NMO. Finally, we should note 
that our calculation is applicable, in principle, to energy loss by any harmonic oscillator 
mode coupled to any two-level system, provided that the model Hamiltonian describing the 
interactions and energy losses is that used in the present work. 

Next, we briefly comment on the initial conditions used in these calculations. We have 
assumed that the NMO begins in a Fock state. In some cases, a more appropriate initial 
condition might be one in which the initial state of the oscillator is a mixture of Fock states 
weighted according to a suitable temperature. If this calculation were carried out, it might 
represent a study of the rate of cooling of an NMO, initially at some temperature T, if it 
were placed in a thermal bath at a much lower temperature and also coupled to an NMO 
at a similarly lower temperature. Another possible initial state might be a coherent state of 
the oscillator, i. e., an eigenstate of a. Both initial states may be achievable experimentally. 
Calculations starting from such mixed states would be more demanding numerically, using 
the present approach, but will probably lead to qualitatively similar results regarding the 
rate of energy loss of both the NMO and the CPB. 

Finally, we comment on our use of the RWA. The RWA is beheved justified when the 
coupling between the NMO and CPB is reasonably small, compared to Itjjjq, and also provided 
the NMO and CPB are in resonance. Both conditions are satisfied in the present calculation. 
We have carried out similar RWA calculations off resonance, with results similar to those 
shown in Figs. 1-3, but we believe that the RWA approximation is less accurate in this 
case. To go beyond the RWA leads to a considerably more comphcated set of stochastic 
differential equations. To confirm our results, we have, in fact, gone beyond the RWA 
for the parameters used in the calculations discussed above (NMO and CPB in resonance, 
relatively weak coupling), solving eqs. (7) including all the counterrotating terms. In this 
case, it proves more convenient to solve numerically for = exp[—i{H — Hf^ — H^)t/fi\\ip), 
where H is given by eq. (5), then transform back to get 1-0). The resulting time dependent 
averages, for these parameters, are nearly indistinguishable from those shown in Figs. 1-3, 
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showing that the RWA is indeed an excellent approximation, as expected, for this case. 
For substantially stronger coupling between the NMO and the CPB, or for frequencies far 
off resonance, we have not yet succeeded in obtaining a converged solution for {ipit)) when 
the counterrotating terms are included, showing that, once again as expected, the RWA is 
inaccurate in this regime. 

To summarize, we have demonstrated, within the RWA, and using a suitable time- 
dependent Schrodinger equation with stochastic terms, that the rate of equilibration of 
an NMO with a thermal bath can be considerably speeded up by coupling the NMO to 
a CPB. Using the same approach, we have demonstrated that the decay rate of the CPB 
can be easily reduced by coupling it to a high Q factor cavity, such as a nanomechanical 



211 ] using a very different 



oscillator; this latter result is consistent with previous calculations 
approach. Both of these effects may be very useful experimentally. In particular, the use of a 
CPB as a means of removing energy from an NMO may be helpful in observing macroscopic 
quantum effects in such oscillators. 
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fited from the facilities of the Ohio Supercomputer Center. We thank Professor B. R. Trees 
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FIG. 1: (Color online) Calculated a^{t) versus time t, averaged over 100 realizations of the noise, 
for the case of zero coupling to the NMO (g = 0), and for fixed g/(2ir) = -(50 MHz), and three 
different values of the NMO damping parameter k. 
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FIG. 2: (Color online) Expected number of excitations in the NMO, n{t), plotted versus time. n(i) 
is averaged over 100 noise realizations, and is shown for zero coupling (g = 0) and for g/{2Tr) = - 
50 MHz and three different values of the CPB damping parameter 7 as indicated. 
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FIG. 3: (Color online) — ln[n(t)/n(0)] for different initial Fock states of the NMO. In all cases, 
g/{2-K) = -50MHz, 7 = 0.05(ns)~^ and k = 0.0005(ns)-^ 
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